ANALYSIS OF THE OPTIMAL EXERCISE BOUNDARY OF 
AMERICAN OPTIONS FOR JUMP DIFFUSIONS 



ERHAN BAYRAKTAR AND HAO XING 



Abstract. In this paper we show that the optimal exercise boundary / free boundary of the American put option 
pricing problem for jump diffusions is continuously differentiable (except at the maturity). This differentiability 
result has been established by Yang ct al. (European Journal of Applied Mathematics 17(1):95-127, 2006) in the case 
where the condition r > q + X (e^ — 1) v{dz) is satisfied. We extend the result to the case where the condition 
fails using a unified approach that treats both cases simultaneously. We also show that the boundary is infinitely 
differentiable under a regularity assumption on the jump distribution. 

1. Introduction 

Let (f2,:r,P) be a complete probability space hosting a Wiener process W = {Wt]t > 0} and a Poisson random 
measure N on x R with the mean measure Xdt v(dz) (in which is a probability measure on K) independent 
of the Wiener process. Let F = {^t}te[o,T] be the (augmented) natural filtration of W and N . We will consider a 
Markov process § = {S**; i > 0}, which follows the dynamics 

(1.1) dSt^ ^lSt-dt + a{St-,t)St-dWt + St- [ {e/ - 1) N{dt,dz), 
as the stock price process. We will take /i = r — q + A — A^, in which 

(1.2) [ e'-i^{dz) <oo, 



as a standing assumption. Wc impose this condition on ^ so that the discounted stock prices arc martingales. The 
constant r > is the interest rate, g > is the dividend. The volatility a-{S, t) is assumed to be continuously 
differentiable in both S and t. Moreover, there are positive constants 6 and A such that 

(1.3) 0<S< a{S, t) < A, for aU ^, t > 0. 

Wc should note that at the time of a jump the stock price moves from St^ to St-e^ in which Z is a random variable 
whose distribution is given by v. When Z < the stock price jumps down, when Z > the stock price jumps up. 
In the classical Merton jump diffusion model, Z is a Gaussian random variable. 

In this framework, wc will study the American put option pricing problem. The value function of the American 
put option is defined by 

(1.4) V{S,t)^ sup E{e-''^{K - Sr)+\So = S}, 
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in which So,T-t is the set of stopping times (with respect to the filtration F) taking values in [0, T — t]. The value 
function V is the classical solution of a free boundary problem (see Proposition 2.1). The main goal of this paper 
is to analyze the regularity of the free boundary. We will show that the free boundary is except at the maturity 
T, and C°° with an appropriate regularity assumption on the jump distribution v. For notational simplicity we 
will first change variables and transform the value function V into u and its free boundary s into b (sec (2.6)) and 
state our results in terms of u and b. 

While the continuity of the free boundary of the American put option in jump models has been studied exten- 
sively, for example, by Pham [1997], Yang ct al. [2006] and Lambcrton and Mikou [2008], the differentiability of the 
free boundary was left as an open problem. Even when the geometric Brownian motion is the underlying process 
the differentiability is difficult to establish (see the discussion on page 172 of Pcskir [2005]) and has only recently 
been fully analyzed by Chen and Chadam [2006/07]. In the jump diffusion case, Yang et al. [2006] proved that the 
free boundary is continuously differentiable before the maturity when the parameters satisfy 



When the condition (1.5) is violated, the free boundary of the American option for jump diffusions exhibits a 
discontinuity at the maturity (see Theorem 5.3 in Yang et al. [2006] and equation (3.20) in this paper). This 
behavior of the free boundary was also observed by Lcvcndorskii [2004] and Lambcrton and Mikou [2008] in the 
exponential Levy models. The purpose of our paper is to extend the regularity results of the free boundary to the 
case where (1.5) is not satisfied. We will see that the boundary is differentiable even when (1.5) is violated. 

There arc two critical points in showing the differentiability properties without the condition (1.5): 1) to show 
the Holder continuity of the free boundary, 2) to show that dgV{S,t) is strictly larger than when the point 
{S,t) is close to the free boundary in the continuation region. Wc achieve these two results in Theorem 3.1 and 
Corollary 3.1 respectively. Combining these two properties and a generalization of the result in Cannon ct al. 
[1974] (see Lemma 4.1), wc upgrade the regularity of the free boundary from Holder continuity to continuous 
differentiability in Theorem 4.1. Then wc analyze the higher order regularity of the free boundary making use of a 
technique Scliaeffer [1976] used for the free boundary of a one dimensional Stefan problem on a bounded domain. 

In order to show that the free boundary is continuously differentiable, it is essential that the value function 
y(5', t) is the unique classical solution of the free boundary problem and has a continuous second derivative (see 
(4.5)). In the jump diffusion models, this has been shown by Pham [1997] under condition (1.5). This condition 
was removed in Yang ct al. [2006] and also in Bayraktar [2008]. Moreover, continuous differentiability of the free 
boundary requires the continuity of the cross derivatives of the value function. In the Levy models with infinite 
activity jumps, the value function is not expected to be a classical solution in general. Yet in the literature 
different notions of generalized solutions were explored. For example, Pham [1998] showed that the value function 
is a viscosity solution, Achdou [2008] showed that the value function is the solution in the Sobolev sense and 
Lambcrton and Mikou [2008] proved that the value function is the solution in the distribution sense. Moreover, 
the smooth-fit property (see (2.4)) is also necessary in our analysis (see Theorem 4.1 and equation (5.1)). While 
this property may not hold for general pay-off functions (see Peskir [2007]), it has been shown to hold for the put 
option pay-off in Zhang [1997], Pham [1997] and Bayraktar [2008] in the jump diffusion models. The analysis in 
this paper also applies to the pay-off functions which are continuously differentiable, bounded, convex on [0, +oo) 
and equal to zero in [K,+oo). In fact, the singularity at the strike of the put option pay-off is the source of the 
technical difficulties. Therefore, wc will focus on the put option pay-off in this paper and leave the investigation of 
the boundary behavior for general pay-off functions to future work. 



(1.5) 
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Figure 1 . Our results and the relationships among them. 
A ^ B means that statement A is used in the proof of statement B. 



Section 2: Properties of tlie value functions 

-Lemma 2.4 



I Section 4: Free boundaries are contiiraously differentiable 
Lemma 4.1 




Corollary 3.1 

Section 3: Free boundaries are Holder continuous 



Lemma 4.1 

Section 5: Higher order regularity of the free boundaries 



The rest of the paper is organized as follows: In Section 2, after changing variables we will collect several useful 
properties of the function u, which will be crucial in establishing our main results in the next three sections about 
the regularity of its free boundary. In Section 3, we will introduce an auxiliary function and use it to show that 
the the free boundary is Holder continuous. In Section 4, we will prove the continuous differentiability of the 
free boundary. In Section 5, we will upgrade the regularity of the boundary curve and show that it is infinitely 
differentiable under an appropriate regularity assumption on the jump distribution. Finally, in Section 6, we will 
show that the approximation free boundaries, constructed in Bayraktar [2008], have the similar regular properties 
with the original free boundary. Proofs of some auxiliary results are presented in the Appendix. 

Our main results are Theorems 3.1, 4.1 and 5.1. In Figure 1 wc show the logical flow of the paper, i.e. we show 
how several results proved in the paper are related to each other. 



2. Properties of the value function 

The value function V{S,t) of the American put option for jump diffusions solves a free boundary problem with 
the free boundary s{t). In particular. Theorem 4.2 of Yang et al. [2006] and Theorem 3.1 of Bayraktar [2008] state 
the following: 



Proposition 2.1. V{S,t) is the unique classical solution of the following boundary value problem: 

(2.1) f + 

(2.2) V{sit),t)=K-sit), ie[0,T) 



^ + -a{S, tfS^^ + ^JiS^ -{r + X)V + X I V{Se\tMdz) = 0, S > s{t), 



(2.3) 



V{S,T) = {K - S)+, S>s{T). 



Moreover, the smooth fit property is satisfied, i.e. 



(2.4) 



±V{s{t),t) ^ -1, te[0,T). 
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In the region {{S,t) : S < s(t),t G [0,r)}, V{S,t) also satisfies the following ineguality: 

dV 1 d^V dV f 

(2.5) + -a{S, tfS^^ + I^S— -ir + X)V + X V{Se',tMdz) < 0. 

In the following, let us first change the variables to state (2.1)-(2.5) in a more convenient form: 

(2.6) a; = log(S'), u{x,t) ^ V {S,T - t) , b{t) = log (s {T ~ t)) and a{x,t) ^ aiS,t). 
It is clear from the assumptions ofa{S,t) that 

a is continuously differentiable in both variables and 

(2.7) 

there are positive constants S and A such that < S < a{x, t) < A for all {x, t) G M x [0, T]. 

While the first part of (2.7) will be used in (4.3) and Lemma 4.1, the second part, which makes sure that the differ- 
ential operators involved are uniformly parabolic, will be necessary for Lemma 2.3, Corollary 5.1 and Theorem 5.1. 
For the simplicity of the notation, we will omit the variables of a in the sequel. In terms of the new variables 
introduced in (2.6), (2.1) - (2.5) reduce to the uniformly parabolic boundary value problem 

(2.8) Cu^^-^a'^-(^t,-^a^^^ + ir + X)u-xJ^u{x + z,tHdz) = 0, x>b{t), 

(2.9) u{bit),t) = K - e''^'\ te{0,T], 

(2.10) u(x,0) = (A'-e^)+, x>6(0), 

(2.11) ^u{b{t),t) = -e''^'\ 

(2.12) Cu{x,t)>0, x<b{t),te{0,T]. 

Let us define the continuation region C and the stopping region T) as follows 

C = {{x,t)\ b{t) <x < +00, < t < T}, V = {{x,t)\ - oo < X < b{t), 0<t<T}. 

From Proposition 2.1, it is clear that the boundary value problem (2.8) - (2.10) has a unique classical solution 
u{x, t) in C. 

Remark 2.1. The integral term in (2.8) can also be considered as a driving term, then the integro- differential 
equation (2.8) can be viewed as the following parabolic differential equation with a driving term f{x,t) = A J^u{x + 
z, t)v{dz): 

(2.13) Tt-2^^^- - 2^ j + + ^ 
This "point of view will be useful in the proof of some results in later sections. 

In this section, we will study the properties of u in both the continuation and the stopping regions. Let us 
start from the following proposition from Yang et al. [2006]. It shows that the time derivative of u is continuously 
differentiable across the free boundary. 

Proposition 2.2. dtu{x,t) is a continuous function m R x (0,T]. In particular, for any t G (0,T], 

(2.14) lim —u(x,t) = 0. 

Proof. The proof is given in Theorem 5.1 in Yang et al. [2006], which summarized Lemmas 2.8 and 2.11 in the 
same paper and used a special case of Lemma 4.1 in page 239 of Friedman [1976]. □ 

Moreover, we will show in the following that t —f u{x, t) is strictly increasing function in the continuation region. 



ANALYSIS OF THE OPTIMAL EXERCISE BOUNDARY 



5 



Proposition 2.3. 



(2.15) 



— {x,t)>0, {x,t)eC. 



Proof. The inequality (2.15) is proved in Proposition 4.1 in Yang et al. [2006] using the Maximum Principle for 
the integro-difFerential equations, which can be found in Theorem 2.7 in Chapter 2 of Garroni and Mcnaldi [1992]. 
However, it can be proved using the ordinary Maximum Principle for parabolic differential equations (see Corollary 
7.4 in Liebcrman [1996]). We know that w = ^ satisfies the following equation in C, 



Since w = dtu > in R x (0, T), (2.16) implies that Cvw > 0. If there is a point {xo, to) g C such that w{xQ,to) = 
(i.e. w achieves its non-positive minimum at (xq, to) ), it follows from the strong Maximum Principle that w{x, t) = 
inCnKx{0<t< to}- Together with the fact that w{x, t) = in V, we have that w{x, t) in M x {0 < t < to}. 
As a result, from 

rto 

u{xo,to) — u{xo,0) = / w{xo,s)ds = 0, 
Ja 

we obtain u{xo,to) = (K — e^")^. This contradicts with the definition of the free boundary b{t), because b{to) = 
maxjx e M : u{x,to) = [K — 6^^)+} and xo > b{to). □ 

Combining Propositions 2.2 and 2.3 with the Hopf's Lemma for parabolic integro-differential equations (see 
Theorem 2.8 in page 78 of Garroni and Menaldi [1992]), we obtain that the free boundary is strictly decreasing. 

Lemma 2.1. The function t b{t) is strictly decreasing for t G (0, T]. 

Proof. The proof is given in Theorem 5.4 in Yang et al. [2006]. □ 

In order to investigate the regularity of the free boundary in the later sections, we need more properties of u, 
which we will develop in the following three lemmas. Since the results of these lemmas are intuitive but proofs are 
technical, we will list the proofs of these lemmas in the Appendix A.l. 

It is well known that S — > V{S,t) is uniformly Lipschitz in R_|_ and t —^ V{S,t) is uniformly semi-Holder 
continuous in [0,T] (see Pham [1997]). The following lemma shows the same properties also holds for u{x,t), the 
function that we obtained after the change of variables in (2.6). (The globally Lipschitz continuity with respect to 
x is not a priori clear and one needs to check whether dxuix^t) is bounded.) 

Lemma 2.2. Letu(x,t) be the solution of equation (2.8) - (2.10), then we have 



(2.16) 




(2.17) 



(2.18) 
(2.19) 



\u[x,t)-u[y,t)\ < C\x 
\u{x,t) -u{x,s)\ < D\t 



x,y eR,te [0,r], 
X eR,0 <t,s <T, 



where C and D are positive constants independent of x and t. 



Proof. See Appendix A.l. 



□ 



In the rest of this section, wc will investigate the boundncss of dtu{x, t) and its behavior when x +oo. These 
two properties will be useful to show several results in Sections 4 and 5 (see e.g. (4.2), proof of Lemma 4.1 and 
Remark 5.1). Let us first recall the definition of the Holder spaces on page 7 of Ladyzenskaja et al. [1968]. 
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Definition 2.1. Let ^ be a domain in M, Qt = x (0,T). We denote Qt the closure of Qt- For any positive 
nonintegral real number a, H°'°'/^ {Qt) is the Banach space of functions v{x,t) that are continuous in Qt, together 
with continuous derivatives of the form d^d^v for 2r + s < a, and have a finite norm 

2r+s<[a] 

in which 

= maxQ^\v\, 

2r+s=[a] Q-2<2r+s<Q 

(a\ \V(X, t) — Vix' t)\ 

= sup </?<!, 

(a;,i),(x',t) e Qt ' ' 
\x - x'\ < po 

<v>f' = sup _ ^) o</3<i, 

{x,t),{x,t')€QT ' ' 
\t-t'\ < po 

where po is a positive constant. 

On the other hand, (Q) is the Banach space whose elements are continuous functions f{y) on O that have 
continuous derivatives up to order [a] and the following norm finite 

||/||(.)^^||,./||(0) + |,H;|<"-'"1\ 
i<H 

in which 

y,y'&^,\y-y'\<Po ^ ' 

Here d^f is the jth derivative of f . These Holder norms depend on po, but for different po > 0, the corresponding 
Holder norms are equivalent hence their dependence on po will not be noted in the sequel. 

Using the Holder spaces and regularity results for parabolic equations, we have the following result. 

Lemma 2.3. For any e > 0, dtu[x, t) is uniformly bounded on M x [e, T]. 

Proof. See Appendix A.l. □ 

Remark 2.2. (i) In the statement of Lemma 2.3, t = cannot be included, i.e., \m\t—tQ dtu{x,t) is not uniformly 
bounded in x because dtu — ^a'^d^u + [p — ^cr^) dxU — (r + X)u + A Jj^ u{x + z, t)v{dz) and limf_,o d'^u{x, t) is 

not bounded as a result of non- smoothness of the initial value at x ^ logii'. 

In the following, we will use the previous lemma to analyze the behavior of dtu(x,t) as x ^ +oo. 

Lemma 2.4. 

lim dtu{x,t) = 0, te (0,T]. 
Proof. See Appendix A.l. □ 
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Remark 2.3. Given the result in Lemma 2.3, it is clear from the differential equation (2.13) that d'^u is uniformly 
bounded m M x [e,T], since dxU is uniformly bounded (see Lemma 2.2). Combining with semi-Holder continuity of 
u{x, •) in Lemma 2.2, Lemma 3.1 in page 78 of Ladyzenskaja et al. [1968] now tells us thatd^u{x,-) & H^/^{[e,T]). 
Therefore, combining with the smooth fit property and Proposition 2. 2, we have 

ueC^(Rx (o,r]). 

In the following three sections we will use the properties of the value function we have shown in this section to 
investigate the regularity of the free boundary b{t). 

3. The free boundary is Holder continuous 

3.1. An auxiliary function. Before we begin to analyze the regularity of the free boundary, let us introduce the 
following important auxiliary function, which was also used in Lambcrton and Mikou [2008] to prove the continuity 
of the free boundary in an exponential Levy model: 

(3.1) J{x,t) =qe'' -rK + X [ [uix + z,t) + e''+^ - K]iy{dz), x e R, t e [0,T]. 

Jr 

As a result of the assumption (1.2), J < oo. Moreover, J is closely related to the behavior of the value function u 
in the stopping region, since one can check that 



(3.2) Cu{x,t) = -J{x,t), ioT X < b{t), t e {Q,T], 



(3.3) Cg{x) = Cu{x,0) = - 



qe'^-rK + X [ {e''+' ^ K)^ iy{dz) 
Jr 



= —J{x, 0), for X < log if. 



in which g{x) = (A' — e^)^. As we shall see in the rest of this section, the function J{x, 0) is of special importance. 
We rename it as Jo{x), i.e., 

(3.4) Jo{x) ^qe^'-rK + X ( (e^+" - K)'^ v{dz). 

Jr 

Let us analyze the properties of J . 
Lemma 3.1. (i) J{x,t) > —rK, Yv[n^\^_^ J{x,t) — —rK and lima;|_|_oo ./(a;, i) = +oo, 

(ii) J{x, t) e Ci (M X (0, T]) n C (R X [0, T]) , 

(iii) The functions x J{x, t) and t J{x, t) are non- decreasing. If either either q > or 

(3.5) ^ ((Af, +oo)) > 0, for any M > 0; 

then X — > J{x, t) is a strictly increasing function. On the other hand, if 

(3.6) u((0,oo))>0 

(3.7) 5tJ(x,t)>0, x>b{t),te{Q,T]. 

Proof, (i) The first statement follows from u{x + z, i) > (A^ — e^+^)+ > K — e^+^. The two limit statements follow 
from the Bounded Convergence Theorem. 

(ii) The continuity of u{x, t) on R x [0, T] implies that J is continuous on the same region. For the differentiability. 
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since dxU and dtu are uniformly bounded in R x [e,T] for any e > (see Lemmas 2.2 and 2.3), the Boimded 
Convergence Theorem gives us 

d 



J(x,t)=qe- + X / 
(3.8) "^"^ ^« 



--u(a: + z, i) + e''+'= 
ox 



v{dz) < +00, 



d f d 

TT- J(a;, t) ~ X -77:u{x + z, t)v{dz) < +oo. 



These partial derivatives are also continuous in R x [e,T] as a result of Remark 2.3. Then the statement in (ii) 
follows since the choice of e is arbitrary. 

(iii) It is clear that the functions x — > J{x,t) and t —>■ J[x,t) are nondecreasing functions since x — > u{x,t) + 
and t u{x, t) are nondecreasing. 

The condition (3.5) means that the support of the measure v is not bounded from above. As a result we have 
that the set A = {z : x + z G C} has positive measure, i.e., v{A) > for any a; € R. For any z G A wc have that 
dxu{x + z,t) + e^+^ > 0, which is equivalent to dsV{Se^, t) + 1 > 0. The latter follows from the convexity of the 
function V and (2.4). li z ^ A, then clearly dxu{x + z,t) + e^+^ = 0. Using these facts in the first equation in 

(3.8) , we see that 3.5 yields dxJ{x, i) > in R x [0, T]. On the other hand, when q > the condition assumed on 
u can be dropped. 

Moreover, when x > b{t) (3.6) ensures that iy{A) > 0. Then (3.7) follows from Proposition 2.3. □ 

In the rest of the paper, we will assume either (3.5) or q > and (3.6) are satisfied. Indeed, in the two well-known 
examples of jump diffusions, Kou's model and Merton's model (see Cont and Tankov [2004] p. Ill), in which ly is 
the double exponential and normal distribution respectively, condition (3.5) is fulfilled. 

As the consequence of Lemma 3.1, the level curve 

(3.9) B{t) ^ {x : J{x, i) = 0, t e [0, T]} . 
is well defined. B(0), which is the unique solution of the integral equation, 

(3.10) Jo{x) = ge^ - rK + X [ (e^+" - K)^ v[dz) = 0. 
will be crucial in describing the behavior of b{t) close to (see Section 3.2). 

Remark 3.1. When r ~ Q, Lemma 3.1 (i) implies that B{t) = — oo. On the other hand, the proof in the following 
lemma tell us that B{t) > b{t). Therefore b{t) — — oo in this case. We will assume r > in the rest of the paper 
to exclude this trivial case. 

This level curve B{t) will be crucial in analyzing the regularity properties of the free boundaries in the rest of 
this section. Let us analyze its properties first. 

Lemma 3.2. (i) B{t) is non-increasing, 

(ii) B{t) e ci(o,r] nc[o,r], 

(iii) B{t) > h{t) for t G (0,T]. Here b{t) is the free boundary m (2.8) - (2.10). 
Proof, (i) The proof follows from Lemma 3.1 (iii). 

(ii) We have the contimuty of B because J{x,t) is continuous and strictly increasing in x (see Lemma 3.1 (ii) and 

(iii) ). Let us focus on the differentiability in the following. It follows from Lemma 3.1 (ii) that J{x,t) is a 
function in R x (0,r]. Moreover, it follows from (3.7) and B{t) > b{t) (which we will prove in the Step 1 in (iii)) 
that 

5tJ(a:,to)L=B(t„) >0, toe(0,To]. 
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Therefore, the ImpHcit Function Theorem impUcs that there exists an open set U containing to such that 

B{t) e C\U). 

Then the statement in (ii) follows after pasting different neighborhoods for all points t e (0, T] together, 
(iii) The proof consists of two steps: 

Step 1: First we show that B{t) > If these is a io S i^iT] such that -B(to) < ^(^0)1 from the definition of 
B{t) and the fact that x J{x, t) is strictly increasing, we obtain J(x, to) > for all x e {B{to), b{to)). Combining 
with (3.2), we have 

Cu{x,to) < 0, for any x e {B{to),b{tQ)), 

which contradicts with (2.12). 

Step 2: Second, we show that B{t) ^ b{t),t e (0, T]. Since b{t) < \ogK (thanks to Lemma 2.1) and t —> B{t) 
is non-increasing, it is clear that B{t) > b{t) for any t G (0,t*) where t* ~T l\ sup{t G M+ : B{t) = \ogK}. Hence 
we only need to focus on the region where B{t) < logK. If there is a to G {0,T] such that B{to) = 6(to)j we will 
derive a contradiction in the following. 

First, let us define the region Q = {{x, t) \ B{t) < x < \ogK, t G (0, T]}. Because of the result in Step 1, C C. 
Hence u{x, t) satisfies 

£vu{x,t) — X / u{x + z,t)iy{dz), (x,t)Gri. 



Let us define ^ = x - B(t), u{(, t) = u{x, t) and g{i, t) = [K - e«+^(*))+ = g{x). In the region Vl = {(^, t)| < ^ < 
\ogK - B{t),t G (0, T]} we have 



(3.11) 



du 1 9 d^u 



dt 2 

since B{t) G C^(0,T]. On the other hand, 



^, + B'{t)--cJ^ 



du 



+ (r + A){t = A / u{£, + z, t)u{dz). 



(3.12) 



Cv~9 = -e«+^(*)i?'(t) + ia2g?+B(t) + L + B' {t) - ]^aA e«+^(*) + (r + A) (/f - e«+^(* 



K v{dz) 



qe'' ' — rK + A 
Therefore, we obtain from (3.11) and (3.12) that 

(3.13) Cv {u - g) (e, t) = qe«+^(*) - rK + X f \ui^ + z,t) + e«+^(')+^ - k] v{dz) = J(C + B{t),t), 

for (^,t) G f2. Note that J{x,t) > when x > B(t). Therefore (3.13) yields 

(3.14) Cviu-g)>0, i^,t)en. 

On the other hand, from our assumption ^0 == b{to) — B{ta) = 0. Moreover, there clearly exists a ball B <zVL such 
that = 0} = (eo,to); 2) (w-5)(e,i) > {u-g){£.aM) = for all (e,t) G B, since (u-g)(^,t) = {u-g){x,t) > 

when X > B{t) > b{t). Now applying Hopf's Lemma (sec Theorem 17 in page 49 of Friedman [1964]) to u — g in 
B, we obtain 

d 



TjT (u-g) (Co, to) > 0, 



(3.15) 

which contradicts with the smooth fit property at (^o? to), i.e., d^{u — g){S,o, to) ~ c>x{u — g){b{tQ), to) = 0. □ 
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Remark 3.2. In the proof of Lemma 3.2 (Hi), the reason we work with the domain instead of the domain fl 
is that may not satisfy the interior ball condition (see Theorem 17 in page 49 of Friedman [1964]), which is a 
crucial assumption of the H op f Lemma. If one can show B(t) S C^, the interior hall condition automatically holds 
for f2 (see Remark in page 330 of Evans [1998]). However, B(t) G does not follow directly from the Implicit 
Function Theorem, because J(x,t) is not expected to be a function in a neighborhood of the point {b(tQ),tQ), for 
any to, as a result of the discontinuity of d'^u{x,t) across the free boundary b{t) (see the following corollary). 

As a corollary of Lemma 3.2 (iii), d^u{x, i) does not cross the free bomidary continuously. 

Corollary 3.1. 

(3.16) ^u{b{t)+,t)^ \mi ■^u{x,t)> ~e^'~'\ te(0,r]. 

OX'^ xib{t) OX^ 

(This is eguivalent to limg^g^f-^ dgV{S,t) > 0,t E [0,T).) 

Proof. On the one hand, since B(t) > b(t) and x J{x,t) is strictly increasing, we have 

(3.17) J{b{t),t)<0, te(0,T], 

On the other hand, from the continuity of u, (2.11), (2.8) and Proposition 2.2, it follows that 
(3.18) 

= lim £u(x, t) = --cr^ lim -rr^uix, t) - -a^e''^'^ - \ ge^^*) - rK + \ / u(b(t) + z,t) + e''(*)+^ - K vidz) 
xib[t) 2 xib(t) dx^ ^ ' 2 \ Jr^ 1 ' ' 

The inequality (3.16) now follows from combining (3.17) and (3.18). □ 

3.2. The behavior of the free boundary close to maturity. We arc ready to analyze the regularity of the 
free boundaries. The continuity of the free boundaries for differential equations with or without integral terms 
have been studied intensively, see e.g. Friedman [1975], Pham [1997], Yang et al. [2006] and Laniberton and Mikou 
[2008]. For the American option in jump diffusions, Pham [1997] showed the continuity of the free boundary under 
the technical condition 



(3.19) r> q + X (e^ - l)i^{dz). 

In Yang et al. [2006], this condition was removed in the proof of the continuity. Moreover, in their Theorem 5.3, 
they showed that 

f logi^, r> q + XL (e/ - l)iy(dz) 

(3.20) 6(0+)^ lim 6(t)=min{logA',i3(0)}= J ' " ' , 7+ ' i m ' 

t-o+ [ B{0), r < q + X J {e/ - l)iy{dz) 



in which B{0) is the unique solution of (3.10). The same result has been shown for the exponential Levy models 
in Lamberton and Mikou [2008] . 

3.3. Holder continuity of the free boundary. In the following, the function Jo{x) in (3.4) and the Maximum 
Principle will play a crucial role in showing that t — )■ b(t) is Holder continuous. 

Lemma 3.3. Let b(t) be the free boundary in Lemma 2.1. For any e > 0, if there exists ^ > such that for any ti 
and t2 satisfying e < ti < t2 < T and <2 ~ < ^ one has 

(3.21) uibiti),t) ^u{b{ti),ti) <C,{t2-ti)", ti<t<t2, 
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in which < a < 1 and is a constant that does not depend on ti and t2, then there exists S' € {0, S] such that 
(3.22) b{ti) ~ b{t2) < C',{t2 - ti)^ , 0<t2-ti<6', 

in which C'^ is another positive constant that is independent of ti and t2. 

Proof. This proof is motivated by Lemma 5.1 in Friedman and Shen [2002]. For any ti and t2 such that e <ti < 
t2 <T and t2 — ti < 6, let us consider the domain D = {{x,t) : b{t) < x < b{ti),ti < t < t2}. (In what foUows, we 
wiU choose ti and t2 close to each other, i.e. we will find an appropriate 5' such that ^2 — ^i < <5'.) Let D be the 
closure of the domain D. 

In the following, we will show that the function 

2 



(3.23) 



X{x) = < VCe{t2-tl)^ +P{x^b{t,)) 



b{t2) <x< b{ti) 



satisfies xi^) i^ {u — g){x, t) on the domain D for suitably chosen positive constant /3. 

It is clear that x{x) = 0, when x < 5(ii) - ^(t2 -ii)^ = C We also havex(fe(ii)) = Cc{t2~ti)" > u{b{ti),t)- 
g{b{ti)) for ti < t < t2 because of the assumption (3.21). On the other hand, x(6(t)) > = u{b{t),t) — g{b{t)). 
Therefore on the parabolic boundary of the domain D, we have that x^'^~ 9- We will show that this holds for all 
{x^t) € D. To this end; we will compare Cx with C{u — g) using the Maximum Principle. Note that x is carefully 
chosen so that it has a continuous first derivative and a bounded second derivative. These properties of x niakes 
the application of the Maximum Principle for weak solutions (see e.g. Corollary 7.4 in Lieberman [1996]) possible. 

First, for (.T,t) G D let us estimate the integral term: 

{ ^/Cl{t2 -ti)^ + I3{x + z- &(ti))}' v{dz) 
{VC^(t2 -ii)^ +/3z}%(dz) 



(3.24) 



^ I + z)v{dz) = A 

< A 

< 2A 



< 2 A [C,{t2 



[C,{t2-t,r+f3^z^] u{dz) 
ti)"+/3'Af] . 

for a sufficiently large positive constant M independent of ti and t2. To obtain the first inequality, we used x < b{ti) 
for {x,t) G D. The third inequality follows, because J^e^i'{dz) < +oo in (1.2) and z is bounded from below. 
With the estimate (3.24), we can calculate Cx inside the domain D. 

1 



Cx{x) 



(3.25) 



> - 



-a p 



(M-^V2)^ 
r + A 



WX^ + (r + \)x 



l{x>4} - A / x{x + z,t)v{dz) 



pH^^y^} - 2A - tiT + 



in which E ^ iB^^Jp: 



> -F{t2-ti)'', 

+ cr'^ + 2AM and F = 2ACe are positive constants. 



Recall that for any e > 0, 6(e) < minjlogivT, i3(0)} and that the strictly increasing function Jq defined in (3.4) 
satisfies Jo{x) < for a; < -8(0). Using these observations and (3.3) it can be seen that for any a; < b{e) we have 



(3.26) 

Now choosing 
(3.27) 



Cg{x) = -Jo{x) > -Jo(&(e)) > 0. 
c = -Jo(&(e)) >0 
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and 5' = min{{^Y^°' , 6} and f3 < y/^, wc have that 

Cx{x){x) > -c > C{u-g){x,t), {x,t) G D. 

Considering 'I' = x — u + g, we have C'^ > in D and 'I' > on the parabohc boundary of D. It foUows from the 
Maximum Principle for weak solutions that > in D, i.e., 

(3.28) x{x)>{u-g){x,t), {x,t) e D. 

Observe that {u — g){x,t) = if a; < ^. For any {x,t) € D, since [u — g){x,t) > 0, we can see that x > ^. This 
gives us 

(3.29) inf 6(t) > &(ti) - ^^(^2 -ti)^, 0<t2-h<S'. 

ti<t<t2 p 

We have shown the free boundary h{t) is continuous and strictly decreasing in Lemma 2.1. Along with this fact, 
the inequality (3.29) gives us (3.22) with = y/Cl/p. □ 

Now we are ready to state the main result of this section. 

Theorem 3.1. Let b{t) be the free boundary in problem (2.8) - (2.10), then for any e > if e < ti < t2 < T , and 
t2 — ti is sufficiently small, then 

(3.30) 6(ti)-6(i2)<Ce(t2-ii)i, , 
in which is a positive constant independent of ti and t2. 

Proof. The proof will follow by applying Lemma 3.3 twice. The first application will show that b{t) is Holder 
continuous with exponent ^. Applying Lemma 3.3 for the second time we will upgrade the Holder exponent to |. 
As a result of Propositions 2.2 and 2.3 for any e > 0, ti and t2 satisfying e < ti < t2 < T wc have that 

du 

(3.31) H(6(ti),t)-H(6(ti),ti) < max — (6(ti),s)(t-ti) <Ci(t2-ii), 

ti<s<t Ot 

where Ci = maX(:<s<T dtu(b(ti), s) is a positive constant. Now as a result of Lemma 3.3, we know that there exists 
a sufficiently small constant 6i G (0, T — e] such that 

(3.32) 6(ii)-^fe) < C;(i2-ti)i 0<t2-ii<5i, 

in which C[ is a positive constant that does not depend on ti, t2 and di. 

It follows from Lemmas 2.8 and 2.11 in Yang et al. [2006] and the Sobolev Embedding Theorem (see also (A-27) 
in Appendix A. 3) that for any a < b < log if and t € [^1,^2]; 



(3.33) 



du du 
^(.,t)-^(a.,t) 



< C* |a; — a;| ^ , x, x G (a, 6), 



in which C is a positive constant that does not depend on t. Taking x = b{ti) and x = b{t) in (3.33) and using 
Proposition 2.2, we obtain 

(3.34) 0<^{b{ti),t)<C\b{t,)^b{t)\-- <C\b{t,)^b{t2)\^, t,<t<t2, 

ot 

where the third inequality follows from b(t) being strictly decreasing in Lemma 2.1. Combining (3.32) and (3.34), 
we get 

(3.35) < -^(6(ti),0 < C2(t2 -ti)^ h<t<t2,0<t2~ti<Si. 

ot 
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As a result 

(3.36) uibiti),t)-u{b{ti),h) < max -^ib{h), s){h - h) < C2{t2 - h)i . 

tl<S<t2 Ot 

Applying Lemma 3.3 for the second time, we know that there exists 62 € (0, Si] such that 

(3.37) b{h) - b{t2) < C,{t2 - h)i , 0<t2-h<d2, 

where is a positive constant that does not depend on ti, t2 and 62- □ 

4. The free boundary is continuously differentiable 

In this section, we will investigate the continuous differentiability of the free boundary. In Theorem 5.6 in Yang 
et al. [2006], the authors have shown that b{t) € C^(0, T], with the extra condition 

(4.1) r>q + X I (e^^- l)i/(dz). 

Thanks to Corollary 3.1 and Theorem 3.1, we can show the continuous differentiability of the free boundary without 
imposing this extra condition. 

Remark 4.1. // condition (4-1) is not satisfied, we can see from (3.20) that there is a gap between limj_to+ b{t) 
and 6(0) = logK . Therefore it is impossible to have b{t) to be even continuous at t ^ 0. But we shall see that it is 
continuously differentiable for all t £ (0, T] . 

Let us consider the time derivative dtu{x,t). Recall that u{x,t) is the solution of (2.8) - (2.10). Using the 
assumption (2.7), the time derivative w = dtu{x,t) satisfies the following partial differential equation 

Cvw = h{x,t), X > b{t), t e iO,T], 
(4 2) «;(6(t),t) = 0, lim w{x,t)^0, ie(0,T], 

w{x,0) = liiTi dtu{x,t), X > 6(0), 

in which 

(4.3) kix, t)-xl l-uix + t).id^) + a . id..) . (£ - I 

When x < b{t), we also have w{x, t) = 0. Given u{x, t) and b{t), (4.2) is a parabolic differential equation for w{x, t). 
In this equation, the boundary conditions for w{x, t) along b{t) and at the infinity follow from Proposition 2.2 and 
Lemma 2.4. 

In order to show the differentiability of the free boundary, we need to study the behavior of g'^g^. u at the free 
boundary (by first making sure that the cross derivatives exist in the classical sense), which is carried out in the 
following lemma. 



Lemma 4.1. (i) As a function oft, -^-g^u{b{t)+,t) = lim^|ft(j) -^-gju{x,t) is continuous on (0,T]. 

dxdt 



(ii) Moreover, the function rP „ u{x, t) is continuous for x > b{t), t S (0,T]. 



This lemma is a slight generalization of the result in Cannon ct al. [1974] to the parabolic intcgro-differential 
equation (4.2). Considering the integral term h in (4.2) as the driving term, this lemma follows from using the same 
technique presented in Section 1 of Chapter 8 in Friedman [1964]. We will postpone this proof to the Appendix 
A. 2. We are now ready to state and prove the main theorem of this section. 

Theorem 4.1. Let b{t) he the free boundary in the boundary value problem (2.8) - (2.10), then b{t) £ C^(0,r]. 
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Proof. First, we will show b{t) is diflFerentiable at to e (0, T]. Let us define p = dlu{b{t^)+, to) + e''^^°\ Corollary 3.1 
implies that p > 0. 

For sufficiently small e > 0, it follows from (2.11) that 

1 



e 

Applying the Mean Value Theorem yields 



^u{b{to), to) - §^<b{U -e),to-e) + e'^*-^ - e^t*""') 



0. 



for some y E (0, b{tQ — e) — 6(to)) and t G (0, e). Letting e ^ in (4.4) and using Lemma 4.1 (ii), we obtain 
(4 5) bjto) - bjto - e) _ ^Xbito)+M 



^0 e £Tuib{to)+,to) + e^i^o)' 

which implies that b{t) is differentiable since p > 0. Moreover, from (2.13) and Proposition 2.2, we have 



, , . ^ 2(r + A) /2(a-r-\) \ 2 , 



which is clearly a continuous function of t on i e (0, T], since b{t) is a continuous function and a-{x, t) is continuous 
from our assumption (2.7). Along with Lemma 4.1 (i), we can see from (4.5) that bit) € C^(0,T]. □ 

5. Higher order regularity of the free boundary 

In the previous section, we have proved that the free boundary b{t) is continuously differentiable. In this section, 
we will upgrade their regularity. Throughout this section, for the simplicity of the notation, we will assume that 
(T is a positive constant. In this case, h{x,t) = A ■§f.u{x + z,t)v{dz)^ which is bounded thanks to Lemma 2.3. 
More generally, if cr = (t(x, t), h{x, t) is given in (4.3). If we assume a{x, t) e C°°(R x [0, T]) with all its derivatives 
bounded and S < a < A for some positive constants S and A, the same arguments in this section can still be 
carried through. Because of Lemmas 2.2 and 2.3, we can see from the equation (2.8) that d^u{x, t) is also bounded 
in R X [e,T] for any e > 0. Hence, h{x,t) is also bounded in this general case. 

First, let us derive an identity for b'{t). Since b{t) is differentiable, taking derivative with respect to t on both 
sides of (2.11), we have 

(5.1) + ^ -e''*^^'(*)- 

The term dlu{b{t)+,t) can be represented as 

Plugging (5.2) back into (5.1) and recaUing w = dtu, we obtain 

We can see from equations (4.2) that w{x^t) is the solution of a formal Stefan problem in the unbounded 
continuation regions C. Schaeffer [1976] gave a proof of the infinite differentiability of the free boundary of a one 
dimensional Stefan problem in a bounded domain. By introducing the new variable ^ = he reduced the 
problem into a fixed boundary problem on a bounded domain. However, if we apply the same change of variables 
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we will have unbounded coefficients in the corresponding fixed boundary problem. Instead, similar to the change 
of variables in the proof of Lemma 3.2 (iii), we will define 

£,=x-b{t), v{(„t) =w{x,t), 

in which b{t) is the free boundary in (2.8) - (2.10). The function v{^,t) satisfies the following fixed boundary 
equation, 

(5.4) 57 - 2^"' " + ^'^^ of) a| + + = Me + h(t),t), (e, t) e (0, +oo) X (0, T], 

(5.5) w(0,<)=0, ie(o,r], 

(5.6) z;(e,0)=«;(e + 6(0),0), C > 0. 
Moreover, we have the following identity 

Remark 5.1. Since b{t) e C^{0,T], so for any e > 0, b'[t) is continuous and bounded in [e,T]. On the other 
hand, since dfU is bounded by Lemma 2.3, so /i(^ + b{t),t) = A J^dtu{£^ + b{t) + z,t)v{dz) is also bounded when 
(£,,t) E [0, +oo) X [e,r]. As a result, it follows from Theorem 2.6 in page 19 of Ladyzenskaja et al. [1968] that the 
parabolic differential equation (5.4) with the initial condition v{S.,e) = w{^ + b{e), e) instead of (5.6) has at most one 
bounded classical solution. It follows from the proof of Lemma 4-1 ft) that dtu{x,t) is a bounded classical solution, 
so it is the unique bounded solution of (5.4). 

The following result for parabolic differential equations will be an essential tool in the proof of the main result 
in this section. 

Lemma 5.1. Let us assume w{£,,t) G i7^"'"([0, +oo) x [(5,r]) (for some a and S > 0) satisfies the following 
equation 

(5.8) ^^a^+e^+cw^dJj{^ + z,tHdz), {^,t)em+^)x{5,T)) 

(5.9) w{0,t) = g{t), t e [S,T]. 

We assume thatd /jj (j){£,+z, t)v{dz) e -ff^"^"([0, +oo)x [5, T]) and that coefficients a, £, c also belong to iJ^°'"([0, +oo)x 
[S,T]) with S < a < A for some positive constants S and A, moreover g(t) £ H^^"{[S,T]). Then u)(^,t) G 
jy2+2a,i+a([o,+oo) X [d',T], for any 5' > S. 

Proof Consider a cut-off function T]{t) e Cg°((0, T]), such that r]{t) = when t £ (0, 5] and T]{t) = 1 for t e [6', T]. 
The function w{(,,t) = ri{t)w{(,,t) satisfies 

dw d'^w dw f dn 



w{0,t) = ,^{t)git), t&[S,T], 
From our assumptions we have that 



d I Vitm + ^, tHdz) + ^w{^, t) e H'"-"{[0, +oo) X [S, T]), 
jR cn 

il{t)g{t)eH'+''{[S,T]). 
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Moreover, the coefRcients of the above differential equation are all inside space lf^"'"([0, +00) x [S, T]). In addition, 
this equation is uniformly parabolic as the result oi < S < a < A. It follows from regularity estimation 
for parabolic differential equations (see Theorem 5.2 in page 320 of Ladyzenskaja et al. [1968]) that w{£,,t) G 
ij2+2"'i+"([0,+oo) X [6,T]), which implies w{^,t) e H^+^°''^+°'{[0,+oo) x [S\T]) by the choice ofr]{t). □ 

Remark 5.2. We will apply the previous lemma to w{x,t) ~ dtu{x,t). Because the initial condition for w{x,t), 
\mit-^Qdtu{x,t), is not smooth, we can not apply Theorem 5.2 in page 320 of Ladyzenskaja et al. [1968] to upgrade 
the regularity of w directly. This is the reason we work with w in the proof of the previous lemma. 

In order to apply Lemma 5.1 to (5.4) - (5.7), we need Holder continuous coefficients and value functions. Let us 
first show that the coefficients in equation (5.4) arc Holder continuous. 

Lemma 5.2. Let b(t) he the free boundary m (2.8) - (2.10). Then h(t) e H^+°'{[5,T]) with < a < ^ for any 
6>0. 

Proof. For any (5 > 0, since b{t) E C^{0,T] by Theorem 4.1, the coefficients in equation (5.4) are bounded and 
continuous in [(5, T]. On the other hand, because dtu(x,t) is bounded in M x [S,T] by Lemma 2.3, the function 
h{£. + b{t), t) = A ^u($ + b{t) + z, t)v{dz) is also bounded when (^, t) e [0, +00) x [5, T]. It follows from Theorem 
9.1 in page 341 of Ladyzenskaja et al. [1968] that equation (5.4) has a unique solution v{^, t) G VF^'^([0, M] x [6, T]) 
for any g > 1 and M > 0. 

By the Sobolev Embedding Theorem (sec, for example, Theorem 2.1 in page 61 of Ladyzenskaja et al. [1968]), 
for g > 3, we have v{£„ t) e iJ/'-'^/^Qo, M] x [5, T]) with /3==2-|(l</3<2). asa resuh, we have 

(5.10) ^„(0,i)GH'^([J,r]), withO<^<i. 

Let us analyze the terms in the denominator on the right hand side of (5.7). We have that b{t) G C^{[5,T]) and 
that 

f{b{t),t) = x£u{b{t) + z,tMdz) e C\[6,T]), 
since u{x,t) E C^CRx [6,T]) (sec Remark 2.3). Moreover, this denominator is also bounded away from 0, because 

{fi-r- A)a2e''« + (r + \)K - f{b{t),t) = ^ (^"(KO, t) + e'^'^) > 0, te [5, T], 
where the last inequality follows from Corollary 3.1. It is clear from (5.7) and (5.10) that, 

b'{t)eH'^{[S,T]). 

□ 

As a corollary of Lemmas 5.1 and 5.2, we can improve the regularity of the functions u{x, t). 

Corollary 5.1. Let u{x, t) be the classical solution of the boundary value problem (2.8) - (2.10). Then u{S,+b(t), t) G 
_y-2+2a,i+aQo,+oo) X [5',T]) for any 5' > 0, with a G (0,1/2). 

Proof. Let ^ = x — b{t), K{£,,t) = u{x,t) and (f>{£^ + z,t) = m(^ + b{t) + z,t). Then K{^,t) satisfies a differential 
equation of the form (5.8) and (5.9) in Lemma 5.1 with g{t) = I{ — e'''*) (in fact k satisfies (5.4) when h in the 
driving term is replaced by /). Moreover, by Lemma 5.2, the coefficients in this equation (5.8) are inside space 
H°'{[5,T]) for any 6 > 0, and g{t) G H^+"{[S,T]). In addition, thanks to the assumption (2.7), the equation (5.8) 
is uniformly parabolic. 
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On the other hand, since u{x, t) is uniformly Lipschitz in .x G M and uniformly semi-Holder continuous in t G [0, T] 
(see Lemma 2.2), and h{t) is continuously difFerentiable, it is not hard to see that /jjw(C + ^(i) + z,t)v{dz) G 
ij2"'"([0, +00) X [<5,T]). Moreover, u{^ + b{t),t) G 7J2"'"([0, +00) x [S,T]) again because of Lemma 2.2. Now, the 
statement follows directly from Lemma 5.1. □ 

Armed with Lemmas 5.1, 5.2 and Corollary 5.1, we can state and prove the main theorem of this section. 

Theorem 5.1. Let bit) be the free boundary in (2.8) - (2.10). Assume that v has a density, i.e. v{dz) = p{z)dz. 
Let a G (0,1/2). If p{z) satisfies J^^p{z)dz G H'^"{R-), then b{t) G Hi+°'{[e,T]). On the other hand, if 
p{z) G i?^"i+2"(R_) for i> 1, then b{t) G i?^+5+"([e, T]), for any e > 0. 

Proof. The proof consists of four steps. 

Step 1. From Lemma 5.2 and Corollary 5.1, we have that b{t) G H'^+"{[6,T]) and that u(^ + b{t),t) G 
i/2+2"a+«Qo,+oo)x[5',r])forany(5' > 5> Owitha G (0, 1/2), which implies that (9tu(^+&(t), i) G H'^"'°'{[0,+oo)x 
[S',T]) (see Definition 2.1). 

Step 2. Assume that there is a positive nonintegral real number f3 with 2(3 < 2a + I, such that 



(5.11) 
(5.12) 
(5.13) 



6(t)Gi/i+^([<5,r]), 

ju{i + bit), t) G H^^'^iiO, +CX)) X [S', T]), 
u{^ + b{t),t) G ff2+2/3,i+/3([o^+oo) X [S',T]), 



for S' > 6 > 0. We will upgrade the regularity exponent from /3 to 1/2 + /3, in steps 2 and 3. 

Let us analyze dtu{£^ + b{t),t). For any integers r, s > 0, 2r + s < 2/3, since dtu{£, + b{t) + z,t) = when z < — ^, 
we have 



ga QT 



dC' dt^ Jr dt 



d 



u{^ + bit) + z,t)iy{dz) 



gs gr 



9r dt 



+ 00 



d 



u{S, + b{t) + z,t)p{z)dz 



s-l 



(5.14) 



1=0 

+00 gs gr g 

_^ d^Wdt 



dC'dt^ dt 

+ b{t) + z, t)p{z)dz, 



for any ^ > 0. 

When t is fixed, in the following, we will show 



(5.15) 



gs gr 



d 



dV dt 

For any ^1 > ^2 > such that ~ C2 < PO: we have 
d' 

d^dt' 



i^ + b{t) + z,t)iy{dz) G ij2/3-[2/3]([o,+oo)), for 2r + s ^ [2/3]. 



|-u(ei + b{t) + z, t)v{dz) - — - 
ot 



d 



s-l 



<i{.s>i}E 



(5.16) 



4=0 
+ 00 



gi gr g 



d^" 9r dt 



-6 
-?1 



( "(a + bit) + z,t)~ 7/(6 + b{t) + z, t) 



d^&_d 

gs gr g 

d^Wdt 



u(6 + bit) + z,t)vidz) 

p(-a)-p(-6) 

piz)dz 



u{^i+bit)+z,t) 



piz)dz. 
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Let US analyze the right hand side of (5.16) term by term. When s > 1, since s — 1 < 2f3 — 1 < 2a + ^ — 1, we have 
p{z) e H^^-'^{R_), which implies 



(5.17) 



^{s>l} 



E 



gi gr g 



js-l-i 

rTT-( p(-ei)-p(-6) 



<c\\d, 



'tu 



1(2/3) 



1 2/3- [2/3] 



in which C is a positive constant and || • jj^^'^) is the Holder norm (see Definition 2.1). On the other hand, it follows 
from (5.12) that 

gs gr g 



(5.18) 



+ 00 
-«2 



/ + 00 
P{z)dz<\\d,ur'''\^l~i2\ 
<2 



p{z)dz 

2/3- [2/3] 



Moreover, because p{z) e i/'^-i+2"(]R_) for ^ > 1 or J^^ p{z)dz e i/2"(]R_), we have jZ^p(z)dz G 77^+2"(R_) 
for ^ > 0. In particular, using 2/3 < 2q; + ^, we can see /"^ p{z)dz £ T?-^''"!^'^! (E_). As a result. 



(5.19) 



-?2 
-«1 



<\\dtu\\^'^^ 



■ti(6 + &(0 + z,t) 



p{z)dz 



p{z)dz 



p{z)dz\<C\\dM\^^P'>\il-i2?'-^'^\ 



where C is also a positive constant. Plugging the estimates (5.17) - (5.19) into (5.16), we observe that (5.15) holds. 
When ^ is fixed, using (5.14), it directly follows from (5.11) and (5.12) that 



gs gr 



d 



(5.20) dC^W m""^^ + ^^^^ + ^' ^^""^"^^^ ^ H^~^{[S\ T]), for 2/3 - 2 < 2r + 5 < 2/3. 



Now, (5.15) and (5.20) imply that 



(5.21) 



!-«(? + b{t) + z, t)u{dz) e ff2/3,/3([o^ ^ 



Let u(^,i) be a bounded solution of the boundary value problem (5.4) with the initial condition u(^, (5') 
dtu{S, + b{S'), (5'). The uniqueness in Remark 5.1 implies that 



(5.22) 



v{C, t) - -u{i + bit), t), (e, t) e [0, X [5', T]. 



As a result, the assumption (5.12) implies that 

(5.23) v{^, t) e H^^'l^{[Q, +oo) X [5', T]). 

We wiU apply Lemma 5.1 to (5.4) - (5.6) with ^(C + z,*) = dtu{£, + b{t) + z,t), a = £ = - {p + b'{t) - (j'^/2), 

c = r + X and d = X. Thanks to (5.11), the coefficient I belongs to ( [5, T] ) . The other coefficients already happen 
to reside there since they are constants. Along with (5.21) and (5.23), Lemma 5.1 yields 



(5.24) 

which implies that 
(5.25) 



v{^,t) e i72+2/3.i+,3Qo,+oo) X [S",T]) for any 6" > 5' > 5, 



d 
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and 
(5.26) 



d_ 

di 



+ 6(t), t) e iJ2+2/3,i+/3([o^ +00) X [5", T]), 



by (5.22). 

Using (5.7) and (5.25), wc will improve the regularity of b{t) in the following. From (A-1) we have 



(5.27) 



f{b{t),t)^X / u{b{t) + z,t)i'{dz) 



= A 



-|-oo 



i{b{t) + z, t)iy{dz) + A / {K- e''(*)+")z/(rfz) 



Along with (5.11) and (5.13), we can see from (5.27) that 

(5.28) fm,t)eH^+"{[S",T]). 

Together with (5.11), (5.25) and (5.28), we can see from the identity (5.7) that b'{t) e Hi+f^{[S" ,T]) for any 
5" > 6' . It in turn implies that 

(5.29) 6(0 e i/3+/3([r,T]). 

Step 3. Let us investigate + b{t), t). For any r, s > 0, 2r + s < 2 + 2/3, we have 

QS Qr 



d^' dt 

d" d 



u{^ + b{t) + z,t)v{dz) 



+ 00 QS Q 

m($ + Ki) + t)p{z)dz + 



di" dp- J_ 



u{i + b{t) + z,t)p{z)dz 



Qt Qr 

-^-u[i + bit) + z,t) 



Qi Qr 

^^_^-^^"(e + ^W + ^,t) 

? QS Qr 



j.s-1-i 



di 



for any ^ > 0. It is worth noticing that d^d^u{i + b{t) + z^t)\zi-(^ 7^ 5|9J"u(^ + b{t) + z^t)\z]-(^ for some i and 
Following the same arguments that lead up to (5.21), we can show 



(5.30) 



u{i + b{i) + z, t)v{dz) e if2+2/3,i+/3([o, +00) X [5\ T]), 



given 1 + 2/3 < 2a + £ - 1. 

Now, we can apply Lemma 5.1 to the differential equation u(i + b{t),t) satisfies, taking (5.13) and (5.29) into 
account. This results in 



(5.31) 



L{i + bit), t) e H^+'^'--+^{[0, +CX3) X [S'", T]), 



for any 5"' > 5" . As a result, we have improved the regularities from (5.11), (5.12) and (5.13) to (5.29), (5.26) and 
(5.31), respectively. 

Step 4. For any e > 0, we apply Steps 2 and 3 inductively starting from (3 = am. Step 1. Let n be the number 
of time we apply Steps 2 and 3. Let S'l = S', in which (5' > is as in Step 1. Running Step 2 and 3 once, we obtain 
two constants (5" and 5'(' such that (5.29), (5.31) hold with (3 = a. In the n-th time, n > 2, we choose (5^ — S'J/^i 
and S','/ > S'^ > 5'^, such that 5'^' < e for any 7i so that [e, T] C [C T]. 

The apphcation of Step 2 for the n-th time will give us that b{t) G H^+"+^ {[e,T]). Applying Step 2 for £+1 
and Step 3 for £ times the result follows. □ 
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Remark 5.3. (i) The previous proof has also shown the higher order regularity oj u{x,t), i.e. u{$, + b{t),t) S 

^2+2q+£4+q+2 ^ [e,r]), for any e > 0, under the assumptions of Theorem 5.1, . 

(ii) Note that h{t) e C'^{{Q,T]) without any assumption on the density p{z). If p{z) e i72m-i+2a(]g_-) 
some m>\, then h[t) G iJ 2+™+"([e^ T]). From Definition 2.1 and the arbitrary choice of e, we have that 
b{t) G C"^^^{{Q,T\) under this assumption. 

As a corollary of Theorem 5.1, we have the following sufficient condition for the infinitely differentiability of b{t). 

Corollary 5.2. Let b{t) be the free boundary in (2.8) - (2.10). Assume that v has a density, i.e. v{dz) = p{z)dz. 
If p{z) G C°°(R_) with j^piz) bounded for each £ >1, but not necessarily uniformly, then b(t) G C°°((0, T]). 

Proof. For any m > 1 with p{z) G C^™+^(M_) and derivatives of p{z) up to order 2m + 1 are bounded, it follows 
from Definition 2.1 that p{z) G i?2"-i+2"(R_). As a resuh of Remark 5.3 (ii), we have b{t) G C""+i((0, T]). □ 

Remark 5.4. There are two well-known examples of jump diffusion models in the literature, Kou's model and 
Merton's model (see Cont and Tankov [2004], P-IH): ™ which the density p{z) is double exponential and normal, 
respectively. For both of these densities, it is easy to see that the conditions for Corollary 5.2 are satisfied. Therefore, 
the free boundaries in both models are infinitely difjerentiable. 

6. The boundaries of the approximating free boundary problems introduced by Bayraktar [2008] 

In this section, we want to show that the approximating free boundaries &„(t), constructed in Bayraktar [2008], 
have regularity properties similar to the free boundary b{t). 

Bayraktar [2008] constructed a monotone increasing sequence {u„}„>o that converges to the unique solution 
u{x, t) of the parabolic integro-differential equation (2.8) - (2.10), uniformly. In this sequence, uo{x, t) = {K—e^)^, 
and each Un{x,t) {n > 1) is the unique classical solution of the following parabolic differential equation: 

dUn 1 2<9^'^'?i f 1 2\ ^''^n 



(6.1) CvUn = - -Cr -^-^ - l^fl ^ j + {r + X)Un = fn{x, t), X > bn{t) , 

(6.2) unibnit),t) = K -e''"'^'\ te(0,T], 

(6.3) w„(a;,0) = (i^-e")+, x > 6„(0), 
in which 

(6.4) fn{x,t) = X / Un-i{x + z,t)i^{dz), 



and the free boundary 5„(i) = log (s„ (T — t)) is defined in terms of s,i(-), which is the approximating free boundary 
in Bayraktar [2008]. Moreover, the smooth fit property is also satisfied for each u„, i.e. 

(6.5) ^Wn(&„(t),t) = -e^"(*', tG(0,r]. 
In the region {{x,t)\x < 6„(t), t G (0, T]}, one also has that 

(6.6) CvUn{x, t) - fn{x, t) > 0. 

We can define the approximating continuation regions C„ and the stopping regions T>n as follows 

Cn = {{x,t)\b^{t) <x < +00,0 < t < T}, P„ ^ {{'x,t)\ ~oo <x< &„(t),0 < t < T}, for aU n>l. 

Since {un}n>o is a monotone increasing sequence, the approximating free boundary {bn}n>i is a monotone de- 
creasing sequence. As a result, we have U„>iC„ = C and n„>il?„ = T>. 



ANALYSIS OF THE OPTIMAL EXERCISE BOUNDARY 



21 



The approximating sequences {un}n>i and {6„}„>i have the similar properties with the value function u and 
its free boundary b. Proposition 2.3, Lemmas 2.1, 2.2 and 2.4 have their analogous versions for u„ and 6„ via the 
same proofs only replacing the integral term / in (2.13) by in (6.4). Proposition 2.2 and Lemma 2.3, on the 
other hand, can be slightly modified as follows: 

Proposition 6.1. For all n > 1, 

(i) If dtUn-i{x,t) is bounded m M x [e,T] for any e > 0, then dtUn{x,t) is continuous m M x (0,T] and 

d 

(6.7) lim ^un{x,t) = 0. 

(ii) On the other hand, if lim.j. ^^^(^t^ dtUn{x,t) = for t £ (0,T] and dtUn{x,t) is continuous in W x (0,T], then 
dtUn{x,t) is uniformly bounded m M x [e,T], for any e > 0. 

Proof. Sec Appendix A. 3 for the proof of (i). Under the assumption that limxib„{t) dtUn{x, t) = for t E (0, T], we 
have dtu{x, t) is bounded in the domain {(a;, t) \ 6„(t) < x < Xq, e < t < T} for any e > and ATq > \ogK. Then 
the rest of the proof of (ii) is similar to the proof of Lemma 2.3. □ 

Remark 6.1. To show that assumptions in both (i) and (ii) are satisfied for all Un, n > 1, we need to walk 
through (i) and (ii) successively. Starting from dtUQ{x,t) = (since wo(x,i) = {K — e^)^), (i) tells us that 
lim^lbj(-(') 9tUi(a;, t) = and dtUi{x,t) is continuous m M x (0,T]. Then it follows from (ii) that dtUi{x,t) is 
bounded m M x [e, T] for any e > 0. This result feeds back to (i). Now, as a result of an induction argument it can 
be seen that assumptions in both (i) and (ii) are satisfied for all n. 

Results similar to Lemmas 3.1, 3.2 and Corollary 3.1 can also be shown to hold for each u„, n > 1. Defining 

Jn{x,t) = qe"" - rK + \ I [un-i{x + z,t) + e''+''- - K]v{dz), xeM,te[0,r], 

Bn{t)^{x:Jr,{x,t)^Q,te [0,T]}. 
we obtain the following: 

(6.8) CvUnix^t) ~ X Un-i{x + z,t)v{dz) ^ -Jn{x,t), x < 6„ (t) , t G [0, T] , 

(6.9) X Jn{x, t) is strictly increasing and t Jn{x, t) is non-decreasing for {x, t) gR x [0, T], 

(6.10) B„{t)>bn{t), te{0,T], 

(6.11) lim — M„(a;,t) > -e^-t*), ie(0,r]. 

Moreover, as we can see in the following Proposition, the approximating free boundaries 6„ have the same critical 
value as & at 0. 

Proposition 6.2. For the approximating sequence 6„(t), we have 

f logic, r> q + \L (e^ - l)vidz) 

(6.12) 6„(0+) ^ hm 6„(t) = min{log 5(0)} = J^, ~ \ ^f, . i m ' 

i-o+ 1^ B(0), r <q + X^^^(e' ~l)v{dz) 

in which B{0) the unique solution of (3.10). 

Proof When x < bn{t){t > 0), it follows from (6.6), (6.8) and (6.9) that 

< CDU„{x,t) - A / Un-iix + z,i)y{dz) = -Jn{x,t) < -J„(a;,0) = -Jo(a;). 
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The fact that Jo{B{0)) = and x — *■ Jo{x) is strictly increasing tells us that x < B{0). Hence 6„(t) < B{Q) thanks 
to the choice of x. It is also clear that 6„(i) < \ogK. Then we obtain 

(6.13) 6„(0+) < min{logii:,B(0)}. 

Now, the corollary results from combining (3.20) and (6.13), since {b„}n>i is a decreasing sequence of functions. □ 

Furthermore, the Holder continuity in Theorem 3.1 also holds for 6„, n > 1. In the proof of Lemma 3.3, we only 
need to replace c in (3.27) by min {— 2/(t^ J„(a;, t)\ bn{t) < x < Bn{t), e < t < T} > 0. On the other hand, results 
in Lemma 4.1 also hold for dxtUn, n > 1. Therefore, combining with (6.11), we have from (6.5) that 

Proposition 6.3. 6„(t) e C^{0,T], n > 1. 

Finally, using the following representation 

-fU„{bnit) + ,t) 



ifJ.-r - X)e'''-W + [r + \)K - /„(o„(t), t) 
one can follow the proof of Lemma 5.2 to show that there is a e (0, 1/2) such that 

bnit) e H^+"i[5,T]), for any ,5 >0. 
Appendix A. 

A.l. Proof of Lemmas 2.2, 2.3 and 2.4. 

Proof of Lemma 2.2. The inequality (2.19) is clear, because we have 

\u{x,t)-u{x,s)\ = \V{e'^,T-t)-V{e'=,T-s)\ < D\t ~ s\i . 

In order to prove (2.18), it suffices to check that dxu(x, t) is uniformly bounded in the domain M x [0, T]. Choose 
a constant X > \ogK + 1, we will first prove dxu{x, t) is uniformly bounded in [X, +oo) x [0, T]. Let us consider 
a cut-off function ri{x) G C°°(R), such that 77(2;) = when a; < AT — 1 and 77(2;) = 1 when x> X . Using (2.13) we 
see that v{x,t) = ri{x)u{x,t) satisfies 

C-DV = il{x)f{x, t) + J{x, t), 
v{x,0) =^ ri{x){K ~ 6"")+ , 



where 

(A-1) f{x,t)^x£u{x + z,t)iy{dz), f{x,t)^-^<j'' i^ifu + 2ij'^j - i^^i - ^a' j i/u. 

It is worth noticing that the term rj'dxU in the expression for / vanishes outside a compact domain. Since we also 
have that u{x,t) < K, both f{x,t) and f{x,t) are bounded inRx [0,T]. 

Let G{x,t]y, s) be the Green function corresponding to the differential operator £-p. We can represent v{x,t) 
in terms of G as 

(A-2) v{x,t)= f dyG{x,t;y,0)r^{y){K-ey)+ + f ds f dy G{x,t;y, s) (f{y, s^y) + f{y, s)) . 

Jr Jo Jr ^ ' 

The first term on the right-hand-side of (A-2) will vanish by the choice of ry(?/). On the other hand. Green function 

G{x, t; y, s) satisfies 



\dxG{x,t;y,s)\ < c{t- s) exp -c 
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for some positive constant c, (see Theorem 16.3 in page 413 of Ladyzenskajaet al. [1968]). Since /gdy exp(— c ^^^_^'^ ) < 
d{t — s) ^ for some other positive constant d, we have that 



ds J dy\dxG{x,t;y,s)\< J dsc{t-s)^^^2ct^, 
Using this estimate and the boundness of / and /, the Dominated Convergence Theorem imphes that 

dxv{x,t)^ I ds I dydxG{x,t;y,s){f{y,s)ri{y) + f{y,s)), 



which is uniformly bomidcd. On the other hand, d^v — rfu + rjdxU. By our choice of ?y(a;), we have that dxu{x,t) 
is uniformly bounded on [X, +oo) x [0,T]. 

Moreover, in the stopping region V, we have dxu{x,t) = —e^. This implies that > dxu{x,t) > —e^^*^ > —K. 
On the other hand, since it is continuous dxU is also bounded in the compact closed domain {{x^t)\b{t) < x < 
X, < t < T}. As a result we have that dxu{x, t) is uniformly bounded 

Proof of Lemma 2.3. Let us choose Xq such that Xq > logK. We will first prove that dtu{x,t) is uniformly 
bounded in the domain [Xo,+oo) x [0,T]. Let k{x,t) € C^(M x [0,T]) be such that 

dxk{x,t)\x=Xo = dxu{x,t)\x=Xo, te [0,r], 

and that k{x,0) = 0, x G M.. These two conditions on k are consistent since dxu{x,0)\x=Xo = 0- The function 
v{x, t) = u{x, t) — k[x, t) satisfies 

(A-3) dxv{x,t)\x=x,=Q, 
and 

(A-4) Lvv{x,t) = f{x,t)+g{x,t), x > b{t),t e {0,T], 

in which g{x,t) — —Cx>k{x,t) and / is given by (A-1). Let us define the even extension of v{x,t) with respect to 
the line x = Xq as 



(A-5) v{x,t) 



v{x, t) X > Xq, 

v{2XQ-X,t) X < Xq. 



We similarly define f{x,t) and g{x). From (A-3) and (A-5), we have v{x,t) G C^'^(M x (0,T]) and that it satisfies 
the equation 

Cvv = f{x, t) + g{x, t), {x, i) e M X (0, T], 
v{x,0)=0, xGR. 

Here the initial condition follows from (2.10) and the choice of Xq and k{x, t). 

It follows from (2.18) and (2.19) that f{x,t) is uniformly Lipschitz in x and semi-Holder coninuous in t. So for 
any xi < X2, if we have either X2 < Xq or Xq < xi, then 

fixi,t)- f{x2,t)\<XC{x2~Xi), 

for the same constant C as in (2.18). On the other hand, if xi < Xq < X2, then 

|/(a;i,i)-/(x2,t)| < \f{xi,t)-f{XQ,t)\ + \f{Xo,t)-f{x2,t)\ 

< \C{Xo-xi) + XC{x2-Xq) = XC{x2-xi). 
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As a result of the last two equations we observe that f{x,t) is uniformly Lipschitz in its first variable. It is also 
clear that f{x,t) is semi- Holder continuous in its second variable. Thus, it follows from Definition 2.1 that 

f{x,t) e X [0,T]), for some < a < 1. 

On the other hand, gix,t) € iJ"^"/2(M x [0,T]), because k{x,t) e C^f (M x [0,T]). Combining with the assumption 
(2.7) on cr, the regularity property of parabolic differential equation (see Theorem 5.1 in page 320 of Ladyzenskaja 
et al. [1968]) implies that 

u(x,t) e ij2+"'i+f (R X [o,r]). 

In particular, u{x,t) G iJ2+"^i+"/2([Xo, +oo) x [0,r]). As a result, in [Xo,+oo) x [0,r], dtu{x,t) is uniformly 
bounded by the Holder norm of u(x,t) . Now. the result follows from the continuity of dtuix^t) inside domain 
{(x, t) I h{t) <x<XQ,e<t<T]ioY any e > (sec Proposition 2.2). □ 

Proof of Lemma 2.4. Let > logi^T be the same as in the proof of Lemma 2.3, again choose a cut-off function 
r}(x) G C°°(]R), such that ri{x) = 1 when x > 2Xo and 77(2;) = when x < Xq. Then formally the function 
r]{x)dtu{x,t) satisfies the following Cauchy problem 

Cvw = r]{x)h{x, t) + h[x, i), (a;, t) G R x [to, T], 

where 

h{x,t)^X / dtu{x + z,t)i^{dz), h{x,t) ^ -l-a^ {2ri'd^dtu + r]"dtu) - ( n-l-crA Tj'dtu, 

Jr 2 \ 2 J 

and we choose ri{x)dtu{x,to), for some to g (0, T), as the initial condition. It follows from Theorem 3.1 in page 

346 of Garroni and Menaldi [1992] that this Cauchy problem has an unique classical solution, we call it w. On the 

other hand, we have w{x,t) = r]{x)dtu{x,t). Indeed, it is easy to check that J^^w{x, s)ds is the unique classical 

solution of the Cauchy problem 

Cdv ~ ds (^r]{x)h{x, s) + h{x, s)^ + ri{x)dtu{x,to), v{x,to)^0. 

Note that r]{x) [u{x,t) — u{x,to)] is another classical solution. Therefore w{x,t) = r](x)dtu{x,t) by the uniqueness. 

Using the Green function G{x,t;y, s) corresponding to the differential operator £d, the solution w{x,t) can be 
represented as 

(A-6) w(x,t)= / dyG{x,t;y,to)w{y,to) + ds dyG{x,t;y,s){r]{y)h{y,s) + h{y,s)), 

Jm. Jto JR 

for all {x,t) G M X {to,T]. Since the Green function satisfies 

\G{x, t; y,s)\< C{t - s)"^ exp , (y, s) G M x [0, t). 

The first term in (A-6) is bounded, as long as w{y,tQ) is uniformly bounded. The contribution of rj'dxdtu (in the 
expression for K) to w is given by, 

dyG{x,t;y,s)T]'{y)——u{y,s)^ / dy—[G{x,t;y,s)rf{y)]—u{y,s). 



dyds 7r dy ds 

Now it follows from Lemma 2.3 that both w{x,to) and h{x,t) are uniformly bounded for x £ M.,t G [to,T]. We 
also have that 77' and 77" vanish outside [Xq, 2Xo]. Since lim^^^+oo G{x, t; y,s) = and it can easily be shown that 
limj^^+oo dyG{s, t; y, s) = 0, the Dominated Convergence Theorem implies that 

lim w{x,t)^0, t £ {to,T]. 

X — *+oo 

Then the statement follows from the choice of 77. □ 
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A. 2. Proof of Lemma 4.1. We will first establish a one to one correspondence between solutions of (4.2) and 
solutions of an integral equation of Volterra type. 

Lemma A-1. (i) Let G{x, t; y, s) be the Green Junction associated to the differential operator Cj) and let us consider 
the following nonlinear integral equation of Volterra type, 

(A-7) {l + \<7^m, t)^ v{t) = - J' dsvis)^a^bis),s) 9,G(6(i), t; s) + ^ N,it), to<t<T, 

where Ni{t) J^^^ dy dxG{b{t),t;y,to)w{y,to) and N2{t) = J^^ds J^^^^ dy dxG{b{t),t]y,s)h{y,s). There exists a 
unique solution v to (A-7). The function v{t) is continuous. 

(a) Let w{x,t) be a classical solution of (4-. 2) on [t{),T] with the initial condition w{x,ta) = dtu(x,to), such that 
t dxw{b(t)+,t) is continuous. Then there is a one to one correspondence between w{x,t) and v{t). Moreover 
t) ^ v{t), to<t<T. 

The initial value of (4.2) may not be smooth. This is the reason we take w{x,to) = dtu{x,to), < to < T, as 
the initial condition of (4.2) and consider the differential equation on t G [toi^]. 

Remark A-1. The correspondence in Lemma A-1 is well known for the Stefan problem on heat equation with 
Lipschitz continuous free boundary (see Section 1 Chapter 8 of Friedman [1964]). Along Friedman's line of proof 
we will extend the correspondence to our parabolic differential equation with Holder continuous free boundary. 

Proof of Lemma A-1. Proof of (i). First, because G{b(t),t;b{s), s) and a{b{s),s) are continuous for s G (0,t) (see 
(2.7)), it follows from the classical result on Volterra equations (see Rust [1934]) that the integral equation (A-7) 
has a unique solution v{t) and it is continuous with respect to t G [to, T], as long as Ni(t), i = 1,2, are continuous 
with respect to t. It is not hard to show these functions are indeed continuous, using the continuity of b(t) and the 
following estimates on the Green function G and its derivatives: 

|9fG(x,i;j/,s)| < C(i - s)-^ exp (^-c^^ 
\d:,G{x, t; y, s) - d^G{x, i; y,s)\< C{t - i)'^ {i - s)"^ exp (^-c^— , 
\d^G{x,t-y,s) - diG{i,t-y,s)\ < C|x-5;|"(t-s)-^exp (^-c^^i^— , 

where i = s < t < t, \x" — ?/| = [a; — y| A |i — ?/|, < a < 1, C and c are positive constants. These estimates 
are from Theorem 16.3 in page 413 of Ladyzenskaja ct al. [1968]. 

Proof of (ii) Let us assume that w{x,t) is a classical solution of (4.2). As a result, the following Green's identity 
(see page 27 of Friedman [1964]) is satisfied 

(9/1 9 1 d 

TT n^'^iy^^)G{x,t;y,s)—w[y,s) ~ -cr'^{y,s)w{y,s) — G{x,t;y,s) - w{y,s)G{x,t;y,s)a(jy{y,s) 
r\ o\ "y V2 oy 2 dy 

d d f f 1 \ \ 
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where to<s<t<T,x> b{t) and y > b{s). Integrating both hand side of (A-8) over the domain b{s) < y < +00, 
to < s < t — e, we obtain 



1 . „ . . > /■*"' . 1 



ds hm ^a'^{y,s)dyw{y,s)G{x,t;y,s) - ds -a'^{b{s), s) dyw{b{s)+, s) G{x,t;b{s), s) 

y—^+oc I I 



y 

t — e -1 i*t — e -1 



(A-9) 



ds hm -a {y,s)w(y,s)dyG(x,t]y,s)+ / ds —a (b(s), s) w(b(s), s) dyG(x,t;b(s), s) 

e pt — e 

ds hm w{y,s)G{x,t;y,s)<7ay{y,s)+ / ds w{b{s), s) G{x,t;b{s), s) a(Ty{b{s), s) 

to «^+°° Jto 

+00 

[G'(x, t;y,t- e)w{y, t - e) - G'(x, t; y, to)w{y, to)] 

6(t-c) 
t-e 



ds 

to 



yii™oo ( ^~ ) 'w{y,s)G{x,t;y,s) - [ fJ- - ^cr^(6(s),s) ) w(6(s), s) G(a;, i; 6(s), s) 



ds / dyG{x,t]y,s)h{y,s). 

to Jbis) 

In the seventh term on the left of (A-9), we used wlxjt) = when x < b{t). Using the boundary and initial 
conditions for w{x, t) and the facts that limy^+00 G{x, t; y, s) ~ and lim-y^+00 dyG{x, t; y, s) ~ 0, letting e ^ 0, 
we can write 

w{x,t) = - dsdxw{b{s)+,s) -a^{b{s),s)G{x,t;b{s),s) + / dy G{x,t;y,to)w{y,to) 

Jto ^ Jb{t) 

'^^^^^ + f ds f dyG{x,t;y,s)h{y,s) 

Jto Jb{s) 

= -Mo{x,t) + AIi{x,t) + M2{x,t). 

Before differentiating both sides of (A- 10) with respect to a;, let us recall the jump identity: if p{t), to < t < T, 
is a continuous function and b{t) is the Holder continuous with Holder exponent a > ^, then for every to <t <T, 



(A-11) 



g rt 1 /■* 

^^^s^ / d,s p{s)G{x,t;b{s),s) = -p{t) + / ds p(s) 9^G(x, t; 6(s), s)|^ . 

xlb{t) OX Jf^-^ Z Jj^j 



This identity can be proved in the similar way as in Lemma 1 in Chapter 8 of Friedman [1964]. As commented in 
the paragraph after Lemma 4.5 in Friedman [1975], the proof of Lemma 1 can go through when we replace Lipschitz 
free boundary with Holder continuous free boundary with the Holder exponent a > i. 
Now we will take the derivative of (A- 10) with respect to x to obtain 

(A-12) |.^^,t)^^|_A./,(x,t) 

i=0 ^ 

and let x J, b{t). Since dxw{b{s)+, s) and cr(&(s),s), to < s < t, are continuous and b{t) is Holder continuous with 
exponent a > 5 (see Theorem 3.1), taking p{s) = i(T^(6(s), s) 9j.u'(6(s)+, s) in (A-11), we obtain 

lim ^Mo{t) = lim ^ / dsia2(6(s),s)9,ti;(fe(s)+,s)G(a;,i;6(s),s) 
(A-13) ^ 

= ia2(6(t), t) dMm+, t)+ I ds ].a\b{s),s) dMKs)+, s) dxG{b{t), t; 6(s), s). 

4 Jto ^ 
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On the other hand, by Lemmas 2.2 and 2.3, w{y,to) and h{y,s) are bounded in M x [toj2^]- Using the Dominated 
Convergence Theorem we get 



(A-14) 
(A-15) 



hin^A/i(a;,t) 

xlb{t) OX 

hin^M2(x,t) 

xib{t) ax 



dy d,G{h{t),t- y, to)wiy, to) ^ iVi(i), 



+ 00 

bit) 
t 

ds 

to Jb{s) 



+ 00 



It follows from (A-12) - (A-15) that d^w{h{t)+,t) satisfies (A-7). 

Let us prove the converse. For any solution v{t) of the integral equation (A-7), wc can define w{x, t) as follows 



w{x, t) 



ds v{s) icr^(&(s), s) G{x, t; 6(s), s) + 



+ 00 



dyG{x,t;y,to)w{y,to) 



b(t) 



ds 



to Jb(s) 



+ 00 



(A-16) 



dyG{x,t;y,s)h{y,s), 
to <t <T,x > b{t), 



and w{x,to) := dtu{x,to). We will show in the following that ■w{x,t) is a classical solution of (4.2) and that 
t — > dxw{b{t)+,t) is continuous. 

Now we will show that w{x,t) defined in (A-16) is a classical solution of (4.2) on [^0,2^] with initial condition 
dtu{x, to). By definition w{x, tg) = dtu{x, to). On the other hand we have that limj-^+oo w{x, t) = 0, which follows 
from the facts that lim^r^+oo G{x, t; y, t^) = and a, v{s), w{y, to) and h(y, s) are all bounded. Furthermore, using 
the properties of the Green function and the definition of w (see A-16), we also have that CT>w{x,t) = h{x,t) for 
X > b(t), t G [to^T]. Observe that dtw , dxW and d^w all exist and are all continuous in this domain. 

In the following we will show that dxw{b{t)+,t) = v{t), which implies the continuity of dx'w{b{t)+,t). We 
differentiate w{x, t) with respect to x and let x I b{t). Since v{t) and a arc continuous and b{t) is Holder continuous 
with exponent a > ^, we can apply the jump identity (A-11) with p{s) — ^a^{b{s), s)v{s). Following the steps 
that lead to (A-7) in the first part of the proof, wc obtain 



(A-17) 



dxw{b{t)+, t) = ~-aHb{t),t) v{t) - I dsv{s)\a^{b{s),s) dxG{b{t), t; b{s),s) + ^ N,{t). 



Comparing (A-17) to (A-7), we see that dxw{b{t)+,t) = i'(t); to <t <T. 

Then it remains to show that w{b{t),t) — 0, to <t <T. To this end, since we have already shown C-pw = h, 
w satisfies the Green's identity given by (A-8). Integrating the identity (A-8) and using (A-16) and the fact that 
limj._f4_oo 'w{x, t) = we can write 



(A-18) 



ds w(b{s), s) 



ia2(6(s), s) + aa^(&(s), s)^ dyG{x, t; b{s), s) - (^fi - ^a"^ {b{s) , s)^ G{x, t; b{s),s) 

X > b{t),to <t<T. 



Let X > bit). Integrating both sides of (A-18) on [x, -f-oo) and using the fact that dxG = —dyG, wc obtain 



dsw{b{s), s) 



ds w{b{s), s) 



^a'^{b{s),s) +(T(Tx 



+ 00 



dudxG{u,t;b{s),s) - ( M - -(T'^ib{s),s) 



+ 00 



du G{u, t] 6(s), s) 



-a^{b{s),s) + aax] G{x,t;b{s), s) 



/i- -CT {b{s),s] 



-\-oo 



du G{u^ t; h{s)^ s) 
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Taking the derivative with respect to x, letting x | b{t) and using the jump identity (A-11) with 
p{s) = (^(T^(6(s), s) + aax) w{b{s), s), we arrive at 

(A-19) 

I (^<j\b{s),s)+aax{b{s),syj w{b{t),t) 



to 



^a^{b{s), s) + aa4b{s),s) ] dyG{b{t), t; b{s),s) -(p- ^<j'\b{s), s) ] G{b{t), t; b{s),s) 



dsw{b{s), s) 

Since b{t) is Holder continuous with exponent a > 1/2, we have 



\dyG{b{t),t;b{s),s)\< ^ 



(t-s)*-" 

Therefore both dyG{b{t),t; b{s), s) and G{b{t),t; b{s), s) are integrable. Consequently, it follows from (A-18), (A-19) 
and the Dominated Convergence Theorem that w{b{t), t) = 0, to < t < T. □ 

Proof of Lemma 4-1- Proof of (i). Let v{t) be the unique continuous solution of the Volterra equation (A-7). Define 
w{x,t) as in (A-16). The Lemma A-1 shows that w{x,t) is a classical solution to equation (4.2). Let us define 

u{x,t) = u{x,to) + / w{x,s)ds, X > b{t),to < t < T. 

J to 

It is easy to check that u{x,t) is a classical solution of the equation (2.8) - (2.10) with initial condition u{x,to)- 
Since (2.8) - (2.10) has a unique solution, we conclude that u{x,t) = u{x,t), x > b{t) and to <t <T. Lemma A-1 
also implies that 

dxdtu{b{t)+, t) = dxw{b{t)+, t) = v{t), to<t<T, 
which implies that dxdtu(b(t)+,t), to <t <T, is continuous. The statement follows since to > is arbitrary. 

Proof of (ii). Let {x, t) be such that x > b{t). Choosing to < t such that b{to) < x, we can see that J^^ dsdxG{x, t; 6(s), s) < 
+00. As a result, we have 

d /■* 1 

— Mo(x,t)= j ds-a\bis),s)dxwibis)+,s)dxGix,t;b{s),s). 

We have shown in part (i) that dx'w{b{s)+, s) is continuous with respect to s. It is easy to show dxMo{x,t) is 
continuous around a sufficiently small neighborhood of {x,t). One can also show that the functions dxMi{x,t), 
i e {1,2} are also continuous by similar means. Thus, it is clear from (A-12) that dxdtu{x,t) is continuous in this 
small neighborhood around {x,t). Therefore, the part (ii) of Lemma 4.1 follows, because of the arbitrary choice of 
X and t. □ 

A. 3. Proof of Proposition 6.1 (i). We will use the following result in Lemma 4.1 in page 239 of Friedman [1976]: 

Lemma A-2. For any a < b < \ogK, < ti < t2 < T , if both u{x,t) and dtu{x,t) belong to L'^{{ti,t2); L'^{a,b)), 
then u{t) belongs to C{{ti,t2)',L'^{a,b)). 

In this lemma, L'^{{ti,t2); L'^{a,b)) is the class of maps which map t G (ti,<2) to the Hilbert space L^{a,b). 
On the other hand C{{ti,t2); L'^{a, b)) is the class of continuous maps which map t G {ti, t2) to £^(a, h). 

The proof of (6.7) is similar to that of (2.14): First, we will study the penalty problem associated to the free 
boundary problem (6.1) - (6.5). Then, we will list some key estimates for the solution of the penalty problem. And 
finally using Lemma A-2 we will conclude. We will give a sketch of this proof below. 
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Let us consider the following penalty problem 

+ M< - 9e) = fnix,t), XeR,0<t<T, 

(A-20) 

ul^{x,0) = ge{x), x€M., 

in which < e < 1, ge{x) G C°°(M) such that ge{x) — (K — e^')+ when x satisfies \K — e^| > e. We define 
fn{x, t) = * fn{x, t), where Ce is the standard mollifier in x and t (see Evans [1998] Appendix C4 in page 629). As 
a result, we have f^{x, t) £ C°°(R x (0, T)). Moreover, because fn{x, t) is continuous, f^{x, t) uniformly converge to 
fn{x,t) on any compact domains as e — > 0. On the other hand, from our assumption that dtUn~i{x,t) is bounded 
for any e > and is a probability measure on R, we obtain that 

(A-21) dtfnix, t) is bounded in M x [e, T], for any e > 0. 

Thanks to (A-21), it is easy to see that dtf^{x,t) are uniformly bounded for any e > 0. The penalty functions 
(3eix) is a sequence of infinitely differentiable, negative, increasing and concave functions such that /?e(0) = — Cg < 
— (r + X)K — re. The limit of the sequence is 



lim f3^(x) 



£^0 



0, a; > 0, 
-oo, a; < 0. 



It is well known that the penalty problem has a classical solution (see page 1009 of Friedman and Kinderlehrer 
[1974/75]). Moreover, a proof similar to that of the proof of Theorem 2.1 of Yang et al. [2006] shows that u^{x, t) G 
C°°(]Rx (0,T))nL°°(IR X (0,T)). 

On the other hand, u'^{x, t) satisfy the following estimates for any a <h < \ogK , < ti < t2 < T, 

(A-22) t(^) {x,t)dx<C, te[ii,i2], 



dt 

(A-23, yj^ (^)-^.<<c, 

(A-24, £^'(^<)^^,, + ^'(^)\,,„<ix<C, 

in which C is a constant independent of e. These estimates use similar techniques to the ones used in the proofs 
of Lemmas 2.8, 2.10 and 2.11 in Yang et al. [2006], since fn{x,t) satisfies (A-21). (Similar estimates can also be 
found in Friedman and Kinderlehrer [1974/75]). We will give the proof for the inequality (A-24) below. The other 
inequalities can be similarly obtained. 

Proof of inequality (A-24)- Let us consider Wn{x,t) = dtu^{x,t). Since u^{x,t) S C°°(]R x (0,T)), it follows from 
(A-20) that Wn{x,t) satisfies 

(A-25) CvWn+P',iu:,-ge)wn = ^/^(x,t). 

Let r]{x,t) e C^(]R x (0, T)), such that r]{x,t) = 1 for {x,t) g [a,b] x [^1,^2], and r]{x,t) = outside a small 
neighborhood of [a,b] x [ti,t2]. Multiplying both sides of (A-25) by rf'dtWn and integrating over the domain 
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fit = R X (0, t) in which ti < t < t2, we obtain 

2 



= 



V 

+ (r + A) 



fit 



dxds — 

dWn 



1 2 2 '^^^n dWn , , 



, dWn dWn 



dxds 



fit 



V^P'eWn - ge)Wn^^dxds 



dx dt 

2dWn d 



dxds 



dt dt 



fn{x, s)dxds 



I1+I2 + h + I4 + h + h 



where Ij is the j-th term on the left and a = a{x,t) satisfying the assumption (2.7). In the following, we will 
estimate each Ij separately. In deriving these estimates we will make use of the inequality 



Ia^ + ab + Ib^ > 0, 

6 6 



(A-26) 

for any A, B Cz R. In the following estimations, C will represent different constants independent of e. 

den] dwn dwn 
dx dx dt 



2„2^^n 9u;„^^^^ 



a 7] 



fit 



dx^ dt 



2 2 d^Wn 

dx dxdt 



(77]- 



dxds 



fit 



> 



1 

4 

1 

4 

6^ 



2 2 9 f dwn 
fit''0t[^ 



2 2 
1 V 



f dWr. 



V dx 



dxds 



1 



darj dwn dwn , , 
(77]— -^^dxds 



iit 



dx dx dt 
2 



(a;, t)dx- - / crry— - 



dcji] ( dwn 



fit 



dt \ dx 



dxds 



fit 



dar] dwn dwn , , 

arj— TT-dxds 

dx dx dt 



f dwn 
\ dx 



1 


darj 






2 









dx J \ dxdt J 



> 



4 



/ dWr, 



\ dx 



(x, t)dx — C 



dxds 

1 
6 



dw„ 

~dr 



dxds 



dxds 



f dwr, 



fit 



V dt 



dxds. 



The first four equalities follow from integration by part. The first inequality follows from the assumption (2.7) and 



the inequahty (A-26) with A = and B = 

For (i=3, 4, 5), a similar procedure yields 



. darj dwji 



The last inequality follows from estimation (A-23). 



h>-c 



fit 



/ dWn 

\ dt 



dxds, I4 > — C — 



iit 



dw„ 
"dT 



dxds, h > —C — 



fit 



For Ir, we have 



> 



-C 



2dWnd_ 9 

' dt dt^'' - 6 



d_ 

dt' 



fn dxds 



dw„ 

~dr 



dw„ 



dxds 



dxds. 



fit 



dw„ 
"dT 



dxds. 



The first inequality can be obtained using (A-26), whereas to obtain the last inequality, we use the fact that 
dtfnixjt) is uniformly bounded. Combining all these estimates for Ij, we obtain 



f dwr, 



V dt 



da; as H 

4 



77 



dwn 
dx 



[x, t)dx < C. 



This completes the proof of (A-24). □ 
Using a similar proof to that of Lemma 2.2 of Yang et al. [2006], we can show that u^{x, t) is uniformly bounded. 
Thus there is a subsequence that {Wj/"} converges weakly to m„ in L'^{{a,b); L'^{ti,t2)) for any a < b < \ogK, 
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0<ti<t<t2<T (see Appendix D in Evans [1998] for an account of the concept of weak convergence). 
On the other hand, it follows from the estimates in (A-22) - (A-24) that and are uniformly bounded in 
L^{a, 6), -Q^^ and g"^" are uniformly bounded in t2); L^{a, b)). Therefore there exists a further subsequence 

satisfying 

dUn' dUn d'^Un' d'^Un d'^Un' d'^Un 



dt dt ' dxdt dxdt' dt^ dt^ ' 
where derivatives of u„ arc defined in weak sense (see Appendix D in Evans [1998]). Here, the convergences are 

weak convergences. Since ||u|| < liminfj llun*^^ || (see Appendix D in Evans [1998] ) (A-22) - (A-24) imply that 

^ e L^{{h,t,);L'{a,b)), ^ e L\{t,,t2); L\a,b)). 

Then it follows from Lemma A-2 that the derivative dtUn exists and is inside the space C((ti, t2); L^{a, b)). On the 
other hand, for fixed t G [^1,^2], it also follows from (A-22) and (A-24) and the Sobolev Embedding Theorem (see, 
for example, Theorem 4 in page 266 of Evans [1998]) that 



(A-27) 



— (x,0 -—(:., t) 



< C\x- a;|-^/^, x,x G (a, 5), 



in which C is a positive constant that docs not depend on t. We already know that dtUn{-,t) is a continuous map 
with respect to t, therefore (A-27) implies that 

^eC{{a,b)x{h,h)). 

Therefore dtUn G C (R x (0, T]) because the choice of a, b, ti and ^2 arc arbitrary and 9tu„ G C ([log K, +00) x (0, T]) 

since [logi^, -foo) x (0,r] G C„. Moreover, we have 

fill fjii 
(A-28) lin, ^(x,fo)= lim ^(6(to), t) = 0, 

xib(to) Ot t^t- at 

because (bn{to), t) is inside the stopping region for t < to as bn{t) is decreasing. 

□ 
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